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Abstract. Local maxima at characteristic comoving scales have previously been claimed to exist in the density 
, perturbation spectrum at the wavenumber k — 27r/L]^gg, where iLgg ~ 100-200 ft" ^ Mpc (comoving), at low 

redshift {z ^ 0.4) for several classes of tracer objects, at z « 2 among quasars, and ai z ~ 3 among Lyman 
break galaxies. Here, this cosmic standard ruler is sought in the "lOK" initial release of the 2dF QSO Redshift 
■".^ , Survey (2QZ-10K), by estimating the spatial two-point autocorrelation functions ^(r) of the three-dimensional 

• (comoving, spatial) distribution of the TV — 2378 quasars in the most completely observed and "covered" sky 

regions of the catalogue, over the redshift ranges 0.6 < 2; < 1.1 ("low-z"), 1.1 < 2 < 1.6 ("med-z") and 1.6 < 
Q . z < 2.2 ("hi-z"). Because of the selection method of the survey and sparsity of the data, the analysis was done 

' conservatively to avoid non-cosmological artefacts, (i) Avoiding a priori estimates of the length scales of features, 

I local maxima in ^(r) are found in all three different redshift ranges. The requirement that a local maximum be 

. present in all three redshift ranges at a fixed comoving length scale implies strong, purely geometric constraints 

on the local cosmological parameters, in which case the length scale of the local maximum common to the three 
redshift ranges is 2LLgg — 244 ± 17 Mpc. (ii) For a standard cosmological constant FLRW model, the matter 
density and cosmological constant are constrained to Qm ~ 0.25 ± 0.10, f^A = 0.65±0.25 (68% confidence), 
5— ( ' Qm = 0.25 ±0.15, Qa = 0.60±0.35 (95% confidence), respectively, from the 2QZ-10K alone. Independently of the 

type la supernovae data, the zero cosmological constant model {Q.a = 0) is rejected at the 99.7% confidence level, 
(iii) For an effective quintessence (wq) model and zero curvature, wq < —0.5 (68% confidence), wq < —0.35 
(95% confidence) are found, again from the 2QZ-10K alone. In a different analysis of a larger (but less complete) 
subset of the same 2QZ-10K catalogue, Hoyle et al. (2001) found a local maximum in the power spectrum to exist 
for widely differing choices of f2m and f^A, which is difficult to understand for a genuine large scale feature at 
fixed comoving length scale. A resolution of this problem and definitive results should come from the full 2QZ, 
which should clearly provide even more impressive constraints on fine features in density perturbation statistics, 
and on the local cosmological parameters Sim, and wq. 

Key words, cosmology; observations — cosmology: theory — distance scale — quasars: general — large-scale 
structure of Universe — reference systems 



1. Introduction 

Increasing observational evidence is mounting in favour 
of fine structure (local maxima) in the power spectrum 
or correlation function of the spatial distribution of extra- 
galactic objects, at scales ~ 100-200 Mpc. At redshifts 

Send offprint requests to: B. F. Roukema 



z ~ 2 — 3, these features have been found among quasars 
and Lyman break galaxies (Deng et al. 1994; Broadhurst 
& JafFe 1999; Roukema k Mamon 2000, 2001), and among 
low redshift objects, similar features have been found by 
numerous groups (e.g. Broadhurst et al. 1990; Broadhurst 
1999; da Costa 1992; da Costa et al. 1993; Baugh & 
Efstathiou 1993, 1994; Gaztanaga & Baugh 1998; Einasto 
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ct al. 1994, 1997a,b; Deng, Deng & Xia 1996; Tago et al. 
2001; see Guzzo 1999 for a recent review, or Kirilova & 
Chizhov 2000 for a very extensive reference list). 

Because the scale ~ 100-200 Mpc is well into the 
linear regime of density perturbations and well above the 
present-day turnaround scale at which matter has had the 
time to collapse by self-gravity, these features should occur 
at fixed comoving length scales, provided that the correct 
values of the local cosmological parameters fim (the mat- 
ter density parameter), Qa (the dimensionless cosmologi- 
cal constant) and wq (the effective qiiintessence parame- 
ter, e.g. Efstathiou 1999) are used to convert redshifts and 
angular positions to three-dimensional positions in three- 
dimensional (proper distance, comoving) space of the ap- 
propriate curvature (e.g. Roukema 2001). 

At low redshifts (2;;$ 0.1), peculiar velocities and the 
use of incorrect values of the local cosmological parameters 
make claims either for or against these fine features diffi- 
cult. On the other hand, at high redshifts {l^z ^3), while 
peculiar velocities have a much smaller effect and assiimp- 
tions on the local cosmological parameters are more often 
stated clearly than for low redshift analyses, the main diffi- 
culty is the sparsity of observational data in homogeneous 
and large volume catalogues. 

In Roukema & Mamon (2000, 2001), the objective 
prism survey of quasar candidates by lovino. Clowes & 
Shaver (1996) was analysed, using the a priori claim, 
primarily from the low redshift analyses, that a local 
maximum in the power spectrum or the spatial corre- 
lation function exists at /c = 27r/iyLss, where i^LSS ~ 
130±10 Mpc. In a fourier analysis of the initial release 
of data from the 2dF Quasar Survey (hereafter, 2QZ-10K, 
Groom et al. 2001), Hoyle et al. (2001) claim marginal 
evidence for the existence of a "spike" at 65-90 Mpc, 
depending on the values of the local cosmological param- 
eters. 

Spectroscopic followup of a small part of the lovino 
et al. (1996) survey (P. Petitjean, private communication) 
suggests that it suffers from contamination by stars and by 
quasars at redshifts other than those estimated by lovino 
et al. (1996), so it is possible that the features found were 
either underestimated or overestimated in strength and/or 
shifted in length scale with respect to the true signal. 

On the contrary, the 2QZ-10K consists uniquely of 
diffraction grating spectroscopic redshifts, so should be 
free of this contamination and potentially provide a 
cleaner signal. 

However, extreme care is required in the analysis of 
the 2QZ-10K survey, because of: 

(a) the "angular selection function" , i.e. the very 
anisotropic distribution of this initial release within 
the two survey regions (e.g. fig. 1, Groom et al. 2001); 

(b) the two sources of incompleteness in the regions ob- 
served so far: "spectroscopic incompleteness" (frac- 
tion of objects observed but with unsuccessful spectro- 
scopic identification) and "coverage incompleteness" 



(fraction of target objects in a given field not yet spec- 
troscopically observed); and 
(c) the initial selection of the target catalogue using ubjr 
multi-colour photometry using photographic plates. 

In this paper, a spatial correlation function analysis 

of the 2QZ-10K is performed, as in Roukema & Mamon 
(2001), in order to minimise spurious artefacts from fac- 
tors (a)-(c) and enable a "comoving standard ruler" anal- 
ysis of the data. Moreover, in order to avoid the a pri- 
ori assumption of the scale at which feature(s) defining 
the standard ruler should occur, the constraint is relaxed 
to only require consistency of the scales of feature(s) be- 
tween different redshift intervals of the 2QZ-10K, without 
assuming the scale(s) at which these feature(s) occur. 

As will be seen below, consistency of length scales in 
the comoving spatial correlation function between differ- 
ent redshifts implies constraints on the local cosmologi- 
cal parameters independently of cosmic microwave back- 
ground, type la supernovae, or other constraints on these 
parameters. 

The observational catalogue and method of analysing 
it in a way that deals with points (a)-(c) above are briefly 
described in Sect. 2, results are presented in Sect. 3, the 
question of whether or not the comoving local maximum 
detected is a genuine, cosmological feature rather than just 
a noise fluctuation is discussed in Sect. 4 and conclusions 
are presented in Sect. 5. 

An almost (i.e. perturbed) Friedmann-Lemaitre- 
Robertson-Walker cosmological model is assumed here, 
optionally including a quintessence relation. The Hubble 
constant is parametrised as h = ifo/lOOkm s~^ Mpc"-^. 
Gomoving coordinates are used throughout [i.e. 'proper 
distances', eq. (14.2.21), Weinberg (1972), equivalent to 
'conformal time' if c = 1; see Roukema (2001) for a self- 
contained and conceptually simple presentation of how 
to calculate three-dimensional comoving separations for a 
curved space, by considering these as arc-lengths in four- 
dimensional space]. Values of the density parameter, Cl^a, 
and the dimensionless cosmological constant, fix, are in- 
dicated where used. 

For convenience, the reader is reminded here of the 
dependence of proper distance on Cl^a and Oa, i-e. 



d{z) 



da 



^0 Ji/(i+z) a-y/Oiil7o^^^^K"+^^Aa^ 
where a is the scale factor, 

fl-K = f^m -|- Oa — 1 



(1) 



(2) 



is the (dimensionless) curvature of the observational 
Universe and 



Rc = 



1 



^0 v?^ 



(3) 



is its curvature radius. 

Effective quintessence distance-redshift relations in a 
flat universe are also considered, where Oq = 1 — fim and 
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Table 1. Angular subsamples of the 2QZ-10K survey. Listed 
are limits in right ascension (aj) and declination {5i) in co- 
ordinate degrees; right ascension interval size 6 in great circle 

degrees; approximate solid angle duj in square degrees; num- 
ber of objects A'^; and a rough lower bound to the solid angular 
number density n in units of (sq. deg.)^^. Note that the survey 
is not uniformly complete within these boundaries; complete- 
ness is modelled in the analysis by use of the angular positions 
of the observed quasars. This is why the number density esti- 
mates are only lower estimates. 
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flqa Q is substituted for J^a in Eq. (1), where wq 

is considered to be a constant (e.g. Efstathiou 1999). The 
value WQ = — 1 corresponds to a standard cosmological 
constant metric. 

2. Catalogue and analysis 

2.1. The 2QZ 

The 2QZ-10K is described in Groom et al. (2001). This 
release only includes quasars in regions of 85% or greater 
"spectroscopic completeness" . 

In order to minimise noise generated by "coverage 
incompleteness", only those quasars observed in regions 
with at least 80% coverage are included in the present 
analysis. The six highest solid angular density regions 
are chosen from fig. 1 of Groom ct al. (2001) and listed 
in Table 1, providing six independent sub-samples from 
which to obtain estimates of ^(r). 

Together, these completeness criteria imply a mini- 
mum of about 70% total completeness, apart from the 
incompleteness implied by initial selection of target ob- 
jects using photographic uh^r photometry (as opposed to 
spectroscopic observation of all objects, selected in a single 
filter from GGD photometry). 

In Roukema & Mamon (2000, 2001), the high number 
density, objective prism survey near the South Galactic 
Pole (lovino et al. 1996) was analysed. The right ascension 
and declination subsamples of the survey in the Roukema 
k. Mamon (2001) analysis had solid angular number den- 
sities of 9.8 and 8.8 quasars/sq.deg. respectively. 

The regions of the 2QZ so far completed have solid 
angular number densities similar to this (lower bounds 
are listed in Table 1), but since a larger redshift interval is 
covered, the comoving spatial densities are about a factor 
of three lower. The signal-to-noise ratios of a peak (local 
maximum) in ^(r) may not be as high in the 2QZ as in 
the lovino et al. (1996) survey. Note that the two surveys 



are independent in angular and redshift range, apart from 
a slight overlap. 

2.2. Method 

2.2.1. Calculation of correlation functions 

The correlation functions are calculated in three- 
dimensional curved space via ^(r) = {DD — 2DR/n + 
RR/n'^)/{RR/n'^) where DD, DR and RR indicate num- 
bers of data-data, data-random and random-random 
quasar pairs respectively (Landy & Szalay 1993), and 
n = 20 times more random points than data points are 
used. The random catalogues use (i) the exact set of angu- 
lar positions of the catalogue quasars, and (ii) random per- 
mutations ("z scrambles". Sect. Illb in Osmer 1981) of the 
observational set of redshifts. Bin size is Ar = 5h^^ Mpc. 

In each of the approximately equal redshift intervals 
0.6 < ^; < 1.1 ("low-z"), 1.1 < z < 1.6 ("med-z") and 
1.6 < z < 2.2 ("hi-z"), ^(r) is calculated for each of the 
six angular sub-samples defined in Table 1. 

Within each redshift interval, mean values of ^(r) and 
standard errors in the mean of ^(r) across the six sub- 
samples are calculated, since these are independent sub- 
samples. Each function ^ is smoothed by a Gaussian with 
CT = 15 Mpc. Weighting in inverse proportion to the 
number of pairs in each angular sub-sample is applied in 
order to obtain signal in proportion to the respective num- 
bers of pairs. 

2.2.2. Consistency of spikes at all redshifts 

Although iLss « 130 ± 10 Mpc has been claimed by 
several authors to be the scale of a local maximum in 
the power spectrum or correlation function, this is not 
assumed here. 

Spikes arc detected here as the first local maximum 
in ^(r) above a given minimum separation fj^inj where 
'^min = 1'^'-' '^^'^ 200 h^^ Mpc. In order for the maxima 
not to be too local (e.g. a slightly high value in one bin at 
r surrounded by slightly lower values in immediately adja- 
cent bins at r — Ar and r+Ar), a smoothing by a Gaussian 
of width CT = 15 Mpc is applied when searching for a 
local maximum. 

For a local maximum to be considered detected and 
useful for constraining local cosmological parameters, the 
estimates of its position rjjiax should be consistent in all 
three redshift ranges. Quantitatively, this consistency is 
defined by requiring the values of r^y^g^ to agree within 
Ar = 5/i~^Mpc between the low- 2: and the med-2; sub- 
samples, and within Ar = 5h~^Mpc between the med-2; 
and the hi-z sub-samples, assuming Gaussian errors. 

So, for a given value of r^i^^, the consistency hypothe- 
sis implies rejection probabilities P and confidence levels 
1-P. 

The above test is applied without any criteria regard- 
ing the significance of individual maxima, in order to be as 
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parameter-free as possible, though at the risk of including 
noise. 

For the sake of discussion, an estimate of the signal- 
to-noise ratio of a local maximum can be quantified as 
follows. 

For a maximum at rjjjg^xj define ^_ and ^+ as the 
value of ^ at r^^g,^ and the values of ^ at the first minima at 



and r_|_ below and above 



respectively, and take 



the root mean square value of cr(j)(r) for r G [r_,r_|_]. 
Then, 



(r) 



(4) 



is the definition of signal-to-noise used for comments be- 
low. 



2.2.3. Advantages of z-scrambling 

(i) The angular positions are used because the angular 
selection function of the 2QZ is difficult to quantify by 
any other method, as stated above (a), in Sect. 1. 

(ii) The use of 2;-scrambling (Osmer 1981) is adopted for 
obtaining the 'random' rcdshift distributions, as in 
Roukema & Mamon (2001), in order to avoid redshift 
selection effects, as mentioned above in (c) of Sect. 1 
(c.f. Scott 1991). 

The problems of incompleteness (b) are also minimised 
by this method of generating the random distributions. 

Note that incompleteness is more likely to lead to noise 
than to spurious eff'ects: if a "void" traced by quasars has 
several quasars missing, then it will be poorly traced, but 
will not be changed in diameter. 

This method implies conservative results. Since some 
of the real, very large scale (e.g. ^ Ihr^ Gpc), but small 
amplitude correlations which may occur in the radial di- 
rection occur in both the real and the so-called random 
catalogues, these will cancel out along with the selection 
effect correlations that are (deliberately) cancelled, and 
be lost from the final signal. (These very large scale cor- 
relation cannot occur in the tangential direction in the 
2QZ-10K, because of the angular selection function and 
division into six sub-samples by angle.) 

So, this is an advantage in the sense that selection 
effects are cancelled; but can also be seen as a disadvan- 
tage in that the amplitude of real, non-zero values in the 
correlation fimction, in particular on large scales, may be 
underestimated. Since the phenomenon of interest here is 
the use of local feature(s) (first derivatives) in the correla- 
tion function as a standard ruler for obtaining geometrical 
constraints, the advantage largely outweighs the disadvan- 
tage: the method is conservative. 

While absolute estimates of the correlation fimction 
(or power spectrum) clearly have cosmological relevance, 
the optimal methods for correcting for angular and red- 
shift selection effects are likely to be considerably different 
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Fig. 1. Toy model showing why small scale correlations can 
be underestimated by ^-scrambling when correlated structures 
exist on the scale of the survey volume and parallel to the 
and z directions (Sect. 2.2.4). The left-hand panel shows a 
toy model with filaments, parallel to these two directions, and 
each composed of 20 points. The right-hand panel shows a z- 
scramble of the distribution in the left-hand panel. The two 
"favoured" values of « and z fs in the original fila- 
ments (length units arc arbitrary in the two directions) result 
in many ^-scrambled pairs containing both these "favoured" 
values. Hence, many more close pairs exist in the 2-scrambled 
distribution than in the original, correlated distribution. So, 
relative to the ^-scrambled distribution, the original distribu- 
tion is anti-correlated for pair separations up to the scale of 
the filament thickness, leading to underestimates of ^ on these 
small scales. 
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Fig. 2. The low-2; (0.6 < z < 1.1) spatial two-point auto- 
correlation function ^(r), (Groth & Peebles 1977), for separ 
rations r in comoving units, for a flat universe with (f2m = 
0.3, f^A = 0.7). The a{r) = 15/i~^Mpc Gaussian smoothed 
mean (^) and error bars corresponding to the standard er- 
ror in the mean a"(^) are shown by the thick and thin solid 
lines respectively. The small scale r <, 25-50 ft" ^ Mpc correla- 
tion is underestimated (hence the shading) due to the effect 
of ^-scrambling in a survey of narrow beam width containing 
filaments (see Sect. 2.2.4 and Fig. 1 for details). 
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Fig. 3. The med-2; (1.1 < z < 1.6) spatial two-point auto- 
correlation function ^(r), (Groth & Peebles 1977), for sepa- 
rations r in comoving units, for a flat universe with (f2m = 
0.3, = 0.7), as for Fig. 2. 
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Fig. 5. The low-2; (0.6 < z < 1.1) spatial two-point auto- 
correlation function ^(r), (Groth & Peebles 1977), for separ 
rations r in comoving units, for a flat universe with (f2m = 
1.0, = 0.0), as for Fig. 2. 
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Fig. 4. The hi- 2; (1.6 < z < 2.2) spatial two-point auto- 
correlation function ^(r), (Groth & Peebles 1977), for sepa- 
rations r in comoving units, for a flat universe with (fim = 
0.3, Oa = 0.7), as for Fig. 2. 




100 2^00 
r (in h Mpc) 

Fig. 6. The med-a (1.1 < z < 1.6) spatial two-point auto- 
correlation function ^(r), (Groth & Peebles 1977), for separ 
rations r in comoving units, for a flat universe with (f2m = 
1.0, fiA = 0.0), as for Fig. 2. 



from the method adopted here, and are not considered in 
this paper. 

2.2.4. Small r disadvantage of z-scrambling 

The 2;-scrambling technique can have a disadvantage if 



(i) either selection effect correlations or intrinsic correla- 
tions exist on a scale comparable to that of two of the 
dimensions of the survey vohimc, and 

(ii) these correlations occur as filaments which are in some 
cases aligned with the tangential {&) and radial {£) 
directions. 
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Fig. 7. The hi-z (1.6 < z < 2.2) spatial two-point auto- 
correlation function ^(r), (Groth & Peebles 1977), for sepa- 
rations r in comoving units, for a flat universe with (f2m = 
1.0, Oa = 0.0), as for Fig. 2. 

The effect is an overestimate of the numbers of close 
pairs in the random simulations, on the scale of the "thick- 
ness" of the "filaments" , and, hence, an underestimate of 
^ (r) on scales r smaller than the "filament thickness" . 

This is illustrated in a toy model in Fig. 1, with a 
qualitative explanation in the caption. 

An algebraic explanation of the effect can be made by 
writing the original distribution as points fl in filament 
1 and points in filament 2, at approximately Cartesian 
coordinates {6, z) for i = 1, 10: 

fl « (0,±i) 

fi ~ (±^,0) (5) 

where "»" indicates that the filaments are of about 1 
length unit in thickness. 

The ^-scrambled distribution then contains many 
points 

g' « (0,0) (6) 

which are strongly correlated [around the point (0, 0)] on 
a scale up to about 1 unit, even though the original dis- 
tribution only contained about one point in this region. 

This was not a serious problem for the lovino et al. 
(1996) sample, since all the survey dimensions were well 
above - 100 200 /i-^ Mpc. 

Here, however, the 2QZ-10K patches (see fig. 1 of 
Croom et al. 2001) are about 2 degrees in size, correspond- 
ing to - 50 lOO/i-^Mpc in the 0.6 < z < 2.2 rcdshift 
range, depending on the values of the local cosmological 
parameters. It is likely that at least some genuine filaments 
happen to lie in the radial and/or tangential directions, 



and some of the redshift selection effects (Scott 1991) are 
likely to create statistical excesses of quasars at certain 
favoured redshifts, i.e. creating what could be called "se- 
lection effect filaments" in the tangential direction. 

What is the thickness of a filament, as traced by 
quasars? A minimum thickness is that corresponding to 
a few times the scale of the non-linear correlation length, 
ro ~ 5/i~^Mpc, a maximum thickness is about 50% 
of typical void sizes, i.e. 50% of ~ 50-100 Mpc, i.e. 
r< 25-50/1-1 Mpc. 

This is the scale up to which an underestimate of ^ (r) 
is possible. On larger scales, there should be a smooth 
transition to a slight overestimate of ^(r), due to pair con- 
servation. 

Since small scale correlations are subject to redshift 
evolution and the effects of non-linear growth of struc- 
ture, they would require model-dependent interpretations 
if they were to be used for comoving standard ruler tests. 
So, the possible anti-correlation on small scales, due to z- 
scrambling in a survey with small dimensions, only affects 
the part of ^(r) which is most difficult to use, and should 
not be a problem for the present paper. 

3. Results 

3.1. Correlation functions 

Figs 2-7 show the correlation functions ^(r) for the low- 
z, med-2; and hi-2; redshift intervals of the 2QZ-10K, for 
the local metric parameter values (flm = 0.3, f2 a = 0.7) 
favoured from other observations, and for the now dis- 
favoured Einstein-de Sitter metric (fim = 1-0, = 0.0). 

Keeping in mind that the method of calculation of ^(r) 
adopted here, which uses 2:-scr ambling, is conservative and 
is likely to imply cancellation of some real correlations, i.e. 
to underestimate the amplitudes of the correlations, and 
that either real or selection effect structures in the data 
are likely to cause an anti-correlation at small-scales, as 
indicated by the shading at r 50 /i-^ Mpc in Figs 2-7, 
the following are clear from the figures: 

(i) A local maximum at a fixed comoving scale, r ~ 
240/1-1 Mpc, is clearly present in all three redshift 
ranges for — 0.3, Oa = 0.7), at signal-to-noise 
ratios [defined in Eq. (4)] of {S/N) = 2.4, {SjN) = 1.1 
and (S/N) = 1.6 at low-z, med-2; and hi- 2; respectively. 

(ii) A local maximum at a fixed comoving scale, r « 
ISO/i-^Mpc, is present in the two higher redshift 
ranges for (fim = 0.3, Oa = 0.7), at a signal-to-noise 
ratio of (S/N) = 1.4 at med-z and (S/N) = 1.7 at 
hi-z, but absent in the low-z data. 

(iii) A local maximum at a fixed comoving scale, r « 
90 Mpc, is present in the two higher redshift ranges 
for (Qm = 0.3, r^A = 0.7), but absent in the low-z data. 

(iv) The anti-correlation at r < 50 /j-^ Mpc expected, due 
to "aligned filaments" , the survey geometry and z- 
scrambling, is present at these scales, for both choices 
of local cosmological parameters. This implies a risk 
of systematic error if the r w 90 /j-^ Mpc maximum 





Fig. 8. Confidence intervals for rejecting tiie iiypotliesis tiiat 
the first local maximum in the spatial correlation function at 
r > 100 Mpc occurs at the same comoving position in the 
low-2, med-z and hi-z redshift intervals of the 2QZ-10K, for 
various hypotheses on the values of (fim, ^^a)- The uncertainty 
in the estimate of the position of a local maximum is ALj^^gg — 
5/1"^ Mpc. Rejection levels of 68% > 1 - P, 95% > 1 - P, 
99.7% > 1-P and 1-P > 99.7 are indicated in this figure and 
in Figs 9,10, with gray shadings from white to gray as a visual 
aid. In this figure, the two white squares indicate (f2rn,f^A) 
values satisfying 68% > 1 — P, and nearly all of the figure 
lies in the 1 — P > 99.7 domain. Due to stronger signals in 
Figs 9,10, the contours in these figures are better labelled. A 
thin line indicating a flat universe (f2m + f^A = 1) is indicated 
in this and the following figures. 

were to be used for the consistency of comoving scale 
test, since it could be strongly affected by this effect. 
So, the conservative approach of using Vjj^y^ values of 
100, 150 and 200 Mpc will be retained here. 

(v) The low-2: correlation function for {Qm — 1-0, J^a = 
0.0), Fig. 5, clearly shows local maxima resembling 
those in the low- 2: (fim = 0.3, JIa = 0.7) correla- 
tion fimction, but displaced to smaller separations. In 
contrast, the med-z and hi-z correlation functions for 
(fim = 1.0, JIa = 0.0) only show broad, more or less 
flat features. 

(vi) The best sign of a feature at a consistent comoving 
scale for (f2m = 1.0, Qa = 0.0) is the (S/N) « 2 
maximum at r « 160 Mpc in the low-z corre- 
lation function, and the {S/N) « 0.5 maximum at 
r « 150 Mpc in the med-z correlation function. 
Since low (S/N) has not been used as a criterion for 
excluding maxima, these two features could have, in 
principle, provided an argument in favour of this choice 
of metric parameters, according to the method defined 
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Fig. 10. As for Fig. 8, confidence intervals for rejecting the hy- 
pothesis that the first local maximum in the spatial correlation 
function at r > 200 ft" ^ Mpc in all three redshift intervals. 

above, if the feature were also present in the hi-z sam- 
ple. However, it does not appear to be present in the 
hi-z sample, and as is seen below quantitatively, is 
rejected as a maximum present in all three redshift 
ranges. 
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Fig. 11. As for Fig. 8, confidence intervals for rejecting the 
hypothesis of a local correlation function mgiximum at r > 
200/i~^Mpc present in all three 2QZ-10K redshift intervals, 
for a flat, effective quintessence model, for various hypotheses 
on the values of (f2m,WQ). A thin line where wq = 1, i.e. 
equivalent to a cosmological constant model, is indicated in 
the figure. 



3.2. Comoving consistency of the position of a local 
maximum at all redshifts 

Figs 8, 9 and 10 show the' coiificicncc levels for reject- 
ing the hypothesis that a local maximum is present in all 
three redshift intervals, for r^^in = 100/i^^Mpc, rjjjjn = 
150/i~^Mpc and r^^^ = 200/i~^Mpc, i.e. without pre- 
selecting a particular scale. 



For rjjiin = 100 h ^ Mpc and 



150/i-iMpc, 



only one or two solutions exist. The solution in common 

to both, and also consistent with the test for r^^^ = 
200/1-1 Mpc, is that with (l^m ~ 0.35, « 0.55). 
The rjjijji = 200/1"^ Mpc test provides a much larger 
range of consistent solutions, consistent with the flat, 
(^m = 0.3, JIa = 0.7) metric, but favouring a hyperbolic 
metric [commonly called an "open" metric, either "open 
open" or "closed open"; see Luminct & Roukema (1999) 
or Roukema (2000) to see why an "open" universe may be 
either open or closed]. 

The reason for the difTcrcnt sizes of the 68% con- 
tours in these plots are clear from the correlation func- 
tions (Figs 2-7). The maximum at r » 240/1"^ Mpc for 
(riui = 0.3, Oa = 0.7) is clearly present in all three redshift 
ranges (Figs 2-4), but in the low redshift range for these 
metric parameters (Fig. 2), only slight, noisy ripples, much 
smaller than the error bars, are present above 100 Mpc 




100 2^00 
r (in h Mpc) 

Fig. 12. The \ow-z (0.6 < z < 1.1) spatial two-point auto- 
correlation function ^(r), for an hyperbolic universe with 
(Qm = 0.35, Qa = 0.55), as for Fig. 2. 

and 150 Mpc for matching the r f« 180 Mpc max- 
imum in the med-^; and hi-^r samples. 

Fig. 12 shows why a consistent solution is found at 
{flm ~ 0.35, f2A ~ 0.55); a local maximum in the range 
160 < r < ISO/j-^Mpc is present, though it is not sig- 



nificant given the error bars. Both the 



100 and 



^min = 150 tests consider this peak to be consistent with 
the r « 180 Mpc in the med-2; and hi-^; redshift sam- 
ples. It is clear that this is not a significant solution. 

It is also clear that no local maximum is common to 
all three redshift ranges for the Einstein-de Sitter metric 
(fim = 1-0, Oa = 0.0), nor for a low matter density model 
with zero cosmological constant {ilm = 0.3, ft a = 0.0). 

3.3. What are the implied constraints on the metric 
parameters? 

Given that the r - 180 /i'^ Mpc maximum providing the 
solution in Figs 8 and 9 has a low signal-to-noise ratio, 
and that the r ~ 240 /j-^ Mpc clearly has the strongest 
signal-to-noise ratio in all redshift ranges, the most robust 
constraint on the metric parameters is that of Fig. 10. 
This result is 

Om = 0.25 ±0.10 (68% confidence) 

fim = 0.25 ±0.15 (95% confidence) 

Oa = 0.65 ±0.25 (68% confidence) 

f2A = 0.60 ±0.35 (95% confidence). 



(7) 



The length scale of the local maximum is also con- 
strained. The 68% and 95% confidence intervals imply 
nearly identical estimates and uncertainties in the length 
scale of the local maximum, which can be labelled 2LLsg 
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because of its similarity to twice the length scale claimed 
in many other studies at low redshift. The scale is esti- 
mated by taking the mean and standard deviation of the 
values r^ax of the local maximum for the three redshift 
intervals, either with or without weighting by the rejec- 
tion probability, within either the 1 — P = 68% or the 
1 — P = 95% confidence contour. 
This formally yields 



2LlsS = 242.9 ± 16.1 /(-^Mpc (68% confidence) 
2LlsS = 243.0 ± 16.6 /i^^Mpc (95% confidence) 
2LlsS = 244.4 ± 16.5/1-1 Mpc (68% confidence, weighted) 
2LlsS = 244.3 ± 16.7/1-1 Mpc (95% confidence, weighted) 

(8) 



which can be summarised as 
LlsS = 122±9/i~iMpc (95%). 



(9) 



A similar hypothesis test to the above for rmin = 
200/i-iMpc, i.e. using the r « 240/i-iMpc local max- 
imum, but for a flat metric in an effective quintessence 
model has also been performed (Fig. 11). This yields 

= 0.25 ±0.10 (68% confidence) 
Vtm = 0.2 ±0.2 (95% confidence) 

wq < -0.5 (68% confidence) (10) 
WQ < -0.35 (95% confidence). (11) 

The introduction of external constraints would strengthen 
these results, e.g. flm > 0.2 would imply wq < 
—0.5 (95% confidence), but for this study it is preferred 
to deduce constraints independently of other observational 
analyses. 








0.5 



m 

Fig. 13. Values of the difltcrence statistic D [Eq. (16)], rep- 
resenting the average absolute slope of ^(r) in the interval 
200/i-iMpc < r < SOO/i-^Mpc over all three redshift in- 
tervals, which would depend monotonically on fim and Oa if 
fluctuations in this interval were due only to noise and not 
to a genuine cosmological signal. The maximum value of D is 
D = 1. Contours from D = 0.4 to D = 0.9, in intervals of 
AD = 0.05 are indicated. Shading is darker for lower values 
of _D. It is clear that D does not increase monotonically from 
(fim = 1.1, = -0.1) to (^m = 0.0, ^A = 1.0). The excess 
values of D, ie. the excess average absolute differences, with 
respect to the underlying trend, occur in a region peaked at 
{Qm « 0.25, Qa « 0.55). 



4. Discussion 

Is the local maximum found at 2LLgg = 244± 17 Mpc 
in all three redshift ranges for (fim ~ 0.25, ^Ia ~ 0.65) a 
real, cosmological signal or could it just be noise which 
happens to give the signal expected? Answers to this ques- 
tion can be divided according to whether or not and which 
external information is accepted as a valid prior assump- 
tion. 

There are strong observational justifications for ex- 
pecting that Clm ~ 0.3, CIa ~ 0.7, and there are also 
numerous observational analyses in favour of fine features 
in the power spectrum or correlation of density perturba- 
tions as traced by extragalactic objects, in particular, in 
favour of a local maximum at ^lSS ~ ^'^^ ^ Mpc. 
The former are close to being widely accepted by many 
independent groups, but the latter remain controversial. 

So, the reality of the signal can be discussed in the 
context of 

(i) no assumption regarding Qni, ^^Aj -^^LSS' 

(ii) the assumption that fim ~ 0.3, « 0.7, but no as- 
sumption regarding i^LSS' o'" 



(iii) the assumption that pLgg ~ 130 ± 10 /j-^ Mpc, but 
no assumptions regarding f2rn,i^A; or 

(iv) the assumptions that fim = 0.25, Oa = 0.65 and 
2LlsS = 244 ± 17/1-1 Mpc. 

4.1. (i) No assumption regarding Q,^, Li^ss 

Let us make no assumption regarding Onu^^A and -^LSg- 
Consider the dependence of an arbitrary property of 
^(r), represented by some statistic of ^(r), as a function 
of fim and 

The value of ^(r) in any bin depends on the num- 
bers of pairs falling in the bin, i.e. on the numbers of 
"data-data" (DD) pairs, "data-random" (DR) pairs, and 
"random-random" (RR) pairs. 

For a fixed distribution in redshift and angular posi- 
tion, as fim decreases and/or as flA increases, the proper 
distance separation between any given pair of quasars in- 
creases. So, as rim decreases and/or as VIa increases, the 
full distribution of pair separations is stretched over a 
larger interval of proper separations, and the Poisson error 
per bin increases. This stretching is mostly in the radial 





Fig. 14. Values of the Poisson-corrected difference statistic D' 
[Eq. (17)], as for Fig. 13, but where D is approximately cor- 
rected for the monotonia variation as a function of expected 
Poisson noise, yielding D' . It is clear that excess values of D' 
occur near (fim ~ 0.25, f^A ~ 0.55) and pass through a lin- 
ear degeneracy from this point to (Sim ~ 0.1, Oa ~ 0.3). This 
confirms that a non-Poisson signal is present in ^(r) for a pair 
(fim, f^A) lying somewhere in this degeneracy region. 




Fig. 15. Values of the Poisson-corroctod difference statistic 
D', as for Fig. 14, but for the interval 100/i"^Mpc < r < 
350 Mpc. The excess values of D' seen in Fig. 14 near 
(fim « 0.25, Oa 0.55) remain present. Imperfect correc- 
tion for the Poisson error also yields a background which 
is not quite flat, so that high D' values also appear near 
fim « 0.05, JIa ^0.95, but are due to Poisson error and not 
an excess relative to Poisson error. 



direction, since the radial sizes of the densely observed re- 
gions of the 2QZ-10K are about an order of magnitude 
greater than their tangential sizes. 

This can be shown algebraically as follows. Ignoring 
the angular separations, since these are small, the full pair 
distribution is spread over the interval which is approxi- 
mately 



(12) 



where d(Clui, ^A, z-i) are the proper distances from the ob- 
server to the redshift limits Zi. The number of pairs m in 
any given separation bin is 

^ (13) 



(14) 



771 OC . 

d{Q,ui, ^A, Z2) — rf(Om, ^A, zi) 
So, the fractional Poisson error per bin is 

OC y d(Om, f^A, Z2) - d(Oni, f^A, Zi). 

m 

Thus, the Poisson error per bin increases when Cl^n de- 
creases and/or Ha increases [sec cqs (1), (2)]. 

If ^(r) has no fine scale features, i.e. only noise fluc- 
tuations, then the increase in Poisson error should be the 
only systematic effect of changing Orn and/or Q^, and it 
should be a monotonic efi'ect, in the sense that the highest 
amount of noise should occur for the lowest fim and the 
highest Oa calculations. 



The favoured values of Qm and Q\ found above 
[Eq. (7)] are not the lowest and highest (respectively) 
of the domain investigated. This provides one argument 
against the local maximum at 2LLsg = 244 ± 17 Mpc 
being a noise fluctuation. 

An independent argument is to define a difference 
statistic 

Doinra,nA) = ^\^ir + Ar,nni,^A)-^{r,flm,^A)\ 

ri 

(15) 



and its normalised value 



maxnjjj,nA{Do(fim, ^^a)} ' 



(16) 



over a range of pair separations ri < r < r2. The lim- 
its ri,r2 can be chosen in order to answer the question 
raised above regarding the 2LLgg local maximum, i.e. 
ri = 200 Mpc, r2 = 300 Mpc, or can be chosen to 
further remove prior knowledge of 2iyLgg while remaining 
in a domain where noise signals are least likely to occur, 
i.e. ri = 100/1-1 Mpc, ra 350 /i'^ Mpc. 

This statistic does not contain any information on 
whether or not a local maximum is present in this range. 
It can be said to represent the average absolute value of 
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the slope of ^(r), or more informally, how much ^(r) fluc- 
tuates. 

If only noise is present, then D should vary monotoni- 
cally as a function of fim and Oa, since the slopes of ^(r) 
would only be noise-generated. 

On the other hand, if a cosmological fluctuation in ^ (r) 
is present in the data, then whether or not it's a local 
maximum, a local minimum or just the rising or falling 
slope of one or the other, the values of D near the correct 
values of (rim, ^^a) should be higher than for surrounding 
values of (Om, Oa). 

An approximate correction for the monotonic variation 
in D, due to the decreased numbers of pairs per bin and re- 
sulting increased (fractional) Poisson errors for lower f^m 
and higher Qa values, can be made by using dependence 
of the maximum proper distance, d{Vlxxi, f^A, z), across the 
full redshift range, and replacing D by 



difiui, r^A, 2.2) - d{nm, r^A, 0.6) 



(i(0.0, 1.1, 2.2) -d(0.0, 1.1,0.6) 



-1/2 



(17) 



Fig. 13 shows the dependence of D on Om and Qa. 

Instead of D increasing smoothly and monotonically 
from (l^m = 1.1, ^^a = -0.1) to (l^m = 0.0,17a = 1-0), 
as would be the case if fluctuations were only caused by 
noise, it is clear that D does not vary monotonically with 
riyxi and Qa- 

Moreover, the exceptional region is clearly that of the 

excess values of D surroimding (fJjn ~ 0.25, f^A ~ 0.55). 
It is not possible that the region «0.15,f2A~0.8isa 
dip relative to the noise values of D, since any additional 
component to D not caused by noise must cause a positive 
contribution (absolute values are always non- negative). 

The approximate correction for the monotonic depen- 
dence on Poisson noise [Eq. (17)] shows this excess sig- 
nal even more clearly, in D' , shown in Figs 14 and 15. 
Although the Poisson noise dependence is not completely 
removed, since the correction applied is only approximate, 
it is clear that D' has excessively large values in the linear 
region from (Om ~ 0.25, JIa ~ 0.55) to (Om ~ 0.1, SIa w 
0.3), and that this is the case whether or not a relatively 
small region around the 2LLgg = 244 ± 17 Mpc max- 
imum, i.e. 200/1-1 Mpc < r < 300 /i'^ Mpc (Fig. 14), 
or a relatively large region around the maximum, i.e. 
lOO/i-^Mpc < r < 350/1-1 Mpc (pjg_ ^^^^^ -g ^^^^ 

It is difficult to explain these excess values of D' unless 
there is a cosmological signal for some (rim, ^^a) po-ir in 
the degeneracy region from (fim ~ 0.25, ^Ia ~ 0.55) to 
(Om « 0.1, Qa ~ 0.3). 

Suppose a metric value pair outside of this region, e.g. 
(Qra ~ 1.0, f^A ~ 0.0), is valid. Then, either 

(a) there is a fluctuation of cosmological origin in ^ (r) of 
the 2QZ-10K, or 

(b) there is not, i.e. there is only noise. 

Consider each case. 



(a) If (Qm ~ 10, r^A ~ 0.0) and there is a cosmological 
fluctuation in ^{r), then there must be an excess value 
of D and D' near (Jim ~ 1-0, ^a ~ 0.0), since for any 
wrong value of (r^m, ^^a), the "correlated pairs" which 
would have been placed in a single separation bin, say, 
the i^^ separation bin, ri + {i—l)Ar < r < ri+iAr [for 
(Om ~ 1.0, Oa ~ 0.0)], are diluted over several bins 
{n + («'-l)Ar < r < ri-|-i'Ar}i'=j,i_i,i_2,...,i+i,i+2,... 
and the value of D and D' is decreased. However, there 
is not an excess in D and D' at (Om ~ 1-0, f^A ~ 0.0). 
So, case (a) is wrong. 

(b) If (fim ~ 1.0, JIa ~ 0.0) and and there is no cosmo- 
logical fluctuation in ^(r), then the fluctuations in ^(r) 
are Poisson noise fluctuations. This noise must increase 
monotonically as Om decreases and/or as ^a increases 
[Eq. (14)]. Hence, there must not be an excess of D and 
D' near (fim ~ 0.25, Qa ~ 0.55), since this is far from 
the limit (f2m = 0.0, JIa = 1-0). However, there is an 
excess of D and D' near (fim ~ 0.25, Qa ~ 0.55). So, 
case (b) is wrong. 

Since both (a) and (b) are wrong, (fim ~ 1-0, JIa ~ 
0.0) is invalid independently of whether or not the exis- 
tence of a fluctuation of cosmological origin is assumed. 

The same argument holds for other pairs (Omj^^A) 
outside of the degeneracy region from (Om ~ 0.25, f^A ~ 
0.55) to {flm « 0.1, r^A « 0.3). 

The only reasonable interpretation of Figs 13-15 is 
that {flm, ^a) lies in the degeneracy region from (fim ~ 
0.25,17a « 0.55) to (fim « 0.1, Ha « 0.3) and that $(r) 
has at least some fluctuating component which is a co- 
moving, cosmological signal and not Poisson noise. 



4.2. (ii) The assumption that ilm ~ 0.3, Oa 
no assumption regarding Llss 



0.7, but 



Given the prior estimate that fim ~ 0.3, Ha ~ 0.7, Figs. 2, 
3 and 4 show that a local maximum of signal-to-noise ratio 

[Eq. (4)] of S/N « 3, 1 and 1 exists near 2LlsS = 244 ± 
17 /i"! Mpc in the low, medium and high redshift ranges 
respectively. 

Why should a local maximum due to noise occur at 
the same position in all three redshift ranges, to within 
the precision of one bin (Ar = f)h~^ Mpc), if it were not 
physical? 

The interval between successive maxima is at least 
« SO/j-^Mpc. In the low redshift bin the interval is 
higher, since the larger error bars lead to insignificant 
maxima at lower length scales, but let us be conservative 
and adopt 60 Mpc as the scale along which a random 
maximum could occur. If it is assumed that local maxima 
occur about once within this interval but with random 
phase, then the chance of agreement between two redshift 
ranges to within ±5 /i^^ Mpc is then P ~ 10/60, and be- 
tween all three is P = (10/60)^ « 0.03, i.e. the hypothesis 
of uniform random phases for these local maxima is re- 
jected at the 97% confldence level. 
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The assumption that Om ~ 0.3, Oa ~ 0.7 also suggests 
that the local maxima which appear at r « 100 Mpc 
and r « 180 Mpc in the med-z and hi- 2; samples could 
be real. 

The former is likely to be affected by the existence of 
correlated structures in a survey geometry of narrow di- 
mensions and the use of z-scrambling (Sect. 2.2.4), so is 
unlikely to be cosmological in origin. The latter should be 
considered a better candidate for a real local maximum. 
However, given the lack of its detection in the low-z sam- 
ple, it is more prudent to consider the 2QZ-10K data to 
be insufficient to establish its existence. 

4.3. (Hi) The assumption that 

Llss ~ 130 ± 10 Mpc, but no assumptions 
regarding nm,^A 

If the existence of a local maximum near ^lSS ^ ^"^^ ^ 
10/i~^Mpc is assumed, then in the analysis discussed 
in Sect. 3.2 and shown in Figs 8 and 12, a search for 
a consistent local maximum in the (f2m) ^a) plane for 
rjjiiji = 100 Mpc should have led to a significant so- 
lution. However, only a weak solution exists, and Fig. 12 
shows that this is for a local maximum at 160 < r < 
180 Mpc rather than at 130 ± 10 Mpc, with a neg- 
ligible signal-to-noisc ratio. 

At first sight, this is, therefore, the most puzzling result 
of this analysis. If there are genuine fine features signifi- 
cantly detected in ^(r) of the 2QZ-10K, why should the 
most commonly claimed feature at low redshift be absent? 

Two likely answers include: 

(i) the sparsity of the sample, and/or 

(ii) coincidence of genuine correlations with redshift and 
or angular selection effects 

which are both capable of removing underlying matter 
density correlations from the sample. No amount of cor- 
rection can (validly) re-introduce correlations which are 
absent from a catalogue due to one or both of these rea- 
sons. 

Another possibility is: 

(iii) that by use of 2;-scrambling and observational angu- 
lar positions to mimic selection effects as closely as 
possible, the cosmological correlations have also been 
mimicked in the random catalogues and, hence, been 
cancelled out of ^(r). 

Note that these effects arc related to the use of a con- 
servative technique: careful correction for selection effects 
risks cancelling some real, cosmological correlations, but 
implies that any detected signal is very unlikely to be a 
selection effect. 

If one or several of the explanations (i)-(iii) are correct, 
then it should still be expected that if a local maximum is 
detected at a scale where sparsity or selection effects are 
less serious, it is detected at a harmonic of the assumed 
scale. 



This is in(lec;d the case. Without prior constraints on 
(fim, ^^a), a local maximum at the scale 2LLgg = 244 ± 
17 Mpc was detected. 

This supports the claim that a loc;al maximum in ^(r) 
or the power spectrum exists at ^lsS ^ 130± 10 h~^Mpc, 
and possibly provides a more accurate estimate (^lsS = 
122±9 Mpc) than the low redshift estimates, since the 
effects of peculiar velocities are much smaller than for the 
low redshift samples. 

Earlier arguments related to (i) include the suggestion 
by de Lapparent et al. (1991) that the existence of a scale 
of ^LSS 120 Mpc, traced by positive density fluctu- 
ations, is implied by the observed existence of large scale 
structure at very low redshifts (de Lapparent et al. 1986), 
consisting of negative density fluctuations on a scale of 
around 20-50 Mpc. 



4.4. (iv) The assumptions that ilni = 0.25, J^a 
and 2Llss = 244 ± 17 /i^^ Mpc 



0.65 



Even if the values of the metric parameters are assumed 
and a scale is chosen, there is no reason why a Poisson 
noise peak should occur at this particular scale (claimed 
to be of interest in several observational analyses) in all 
three redshift intervals, in excess of Poisson noise peaks 
generated from uncorrclatcd, random simulations, unless 
it is cosmological in origin, provided that the simulations 
mimick any possible peculiar properties of the catalogue. 

So, 200 random simulations, using the z-scrambling 
technique as above, were performed for each redshift in- 
terval, for each angular sub-sample, for (fim = 0.25, flA = 
0.65). 

Similarly to Roukema & Mamon (2001), the proba- 
bility Pi that a local maximum can occur as close to 
2LLgg — 244 ± 17/1^^ Mpc in the simulation as in the 
observational catalogue, with at least the same signal-to- 
noise ratio [Eq. (4)] as the observational local maximum, 
for the i^^ redshift interval, is deflned 



2iLSSl < kobs-2i^LSsl and 



(S/7V)S™ > (S/iV)0bs 



(18) 



where the position and signal-to-noise ratio of the first 



local maximum at r > 

^obs 



200 h ^ Mpc are r^^jg and 



(5/iV)°°^ in the observations and rg^^ and {S/Nf^^ in 
the simulations. The combined probability for the three 
redshift intervals is 



A 23 = P1P2P3 



(19) 



since the observational data sets in the three redshift in- 
tervals are independent. 

The results are Pi = P2 = 0.065, P3 = 0.075, so that 
P123 =3xl0-f 



In other words, given that {Qm — 0.25, JIa = 0.65), 
the probability that random Poisson signals give a local 
maximum as close to 2LLgg = 244 Mpc as the lo- 
cal maximum in the observational data and of at least 
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as strong a signal-to-noise ratio as the observational local 
maximum, in all three redshift intervals, is rejected at the 
99.97% confidence level. 

4.5. Comparison with the analysis of Hoyle et al. 
(2001) 

How do these results compare with the fourier spectrum 
analysis by Hoyle et al. (2001)? 

Given the large redshift range of the 2QZ-10K (80% 
of the quasars lie in the range 0.6 <, z ;< 2.2), it is puz- 
zling that Hoyle et al. (2001) found a local maximum 
to exist for widely differing choices in the local cosmo- 
logical parameters: it appears to be present both for 
{rim = 0.3,12a = 0.7) and for (J^m = 1.0, f^A = 0.0), 
at r ~ 90/i~^Mpc and r w 65/i~^Mpc respectively in 
their analysis. 

Comparison of fig. 8(a) and fig. 8(b) of Hoyle et al. 
(2001), representing the cases (r^m = 1.0, Qa = 0.0) 
and ( flm = 0.3, il\ = 0.7) respectively, seems to show a 
stronger signal in the latter case, but this is still problem- 
atic for a feature which should exist at a fixed comoving 
scale only for the correct values of the local cosmological 
parameters. 

However, differences in the technique of Hoyle et al. 
(2001) relative to our own include: 

(i) Hoyle et al.'s use of a very smooth redshift selection 
function, which probably does not correct for the se- 
lection effects discussed by Scott (1991); 

(ii) Hoyle et al.'s inclusion of angular selection regions 
that have less than 80% "coverage completeness" (two- 
thirds of their sample lies in regions with only 20%- 
80% coverage completeness), with a statistical correc- 
tion for this via their window function; 

(iii) use of a fourier analysis rather than a correlation 
function; 

(iv) use of large, logarithmic bins in k rather than linear 
bins of Ar = 5 Mpc in r. 

(i) The non-removal of redshift selection effects secims 
a reasonable candidate for inducing a convolution of an 
artefact with the real signal, shifting it to the wrong scale, 
and generating a weak signal for the wrong values of the 
metric parameters. 

(ii) Inclusion of low "coverage completeness" regions 
is likely to introduce strong sources of noise. Even though 
this should be statistically corrected by the window func- 
tion, our own results suggest that it may be more useful 
simply to ignore low "coverage completeness" regions. The 
disadvantage is that only about 3000 of the 11000 quasars 
can be used, but the advantage is having less noise to cor- 
rect. 

(iii) A fourier analysis of data with a very complex 
window function is obviously difficult, and probably more 
useful for estimating the overall shape of the power spec- 
trum than for the detection of fine features in it. Although 
in principle, the power spectrum and the correlation func- 
tion are intrinsically related through Fourier transforms, 



conversion from one to another is diffic;ult in practice and 
will be left to other authors, particularly since a large part 
of the power is removed from our estimate of ^(r) by the 
use of 2;-scrambling: this technique is optimised to detect 
fine features, not gross features. 

It should be noted that in practice, the correlation 
function is more commonly calculated on small scales, 
while the power spectrum is calculated on large scales. 
This is probably because the simplest implementation of 
calculation of the correlation function, via direct counting 
of pairs, is more computationally intensive than that of 
fast fourier transforms. 

However, if (fim, ^^a) parameter space is to be ex- 
plored, then using correlation functions is more precise 
than using Fourier transforms. This is because Fourier 
transforms don't exist in curved space. In curved space 
over the range considered, calculation of Fourier trans- 
forms instead of 3-spherical harmonics or 3-hyperbolic 
eigenfunctions is still a good approximation up to a 
few hundred Mpc, since the curvature radius is at least 
2860 Mpc. Moreover, observations indicate that the 
curvature radius is at least about 10,000 Mpc. 

(iv) The use of wide, logarithmic bins may make it 
difficult to detect a fc^ature which (in tlic^ prc^scnt analysis) 
is much less strong in amplitude than that claimed by 
Einasto et al. (1997b). 

Any or all of these could explain the differences in the 
two analyses. 

It should be noted that Groom et al. (2000) performed 
a correlation function analysis of the 2QZ-10K on scales 
below ~ 100 Mpc, avoiding problem (ui) above, mak- 
ing calculations to consider the effects of (ii) , but retaining 
problems (i) and (iv). Problem (i) is likely to be insignif- 
icant on scales below ~ 100ft~^Mpc, so should not be a 
problem for the purposes of the analysis of Groom et al. 
(2000). 

In this paper, scales above 100 /i""'^ Mpc are those of 
most interest, which is why the correction technique of z- 
scrambling has been used here to avoid problem (i). So, 
the two studies are complementary and cannot be directly 
compared. 

4.6. Comparison with the analyses of Roukema & 
Mamon (2000, 2001) 

As mentioned above (P. Petitjean, private communica- 
tion), some fraction of the lovino et al. (1996) sample 
of high-quality quasar candidates appears to consist of 
quasars at wrongly estimated redshifts and stars. This 
implies that some "redshifts" in the analyses either are 
those of genuine quasars, but for a misidentifed emission 
line, or correspond to the emission wavelengths of atomic 
transitions which should vary little from star to star. 

Both cases should add random noise in the cross- 
correlations, and some genuine correlations should exist 
among the misidentified quasars, but at underestimated 
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tangential separations (if the true redshifts are lower than 
the estimated ones, which is apparently the case). 

However, since the emission lines of the stars should 
be at essentially fixed "false redshifts" (since the velocities 
of stars are much less than cosmological expansion veloc- 
ities), the strongest artefact in the sample is likely to be 
objects at certain favoured "false redshifts" . 

In this case, the purely tangential analysis (Roukema 
& Mamon 2000) is likely to be less affected by contaminat- 
ing objects than the three-dimensional analysis (Roukema 
& Mamon 2001). Since the scale found in Roukema & 
Mamon (2000) is consistent with the present value of 
-^LSS ~ ± 9 Mpc for values of fim consistent with 
those estimated here, the results are compatible. 

4. 7. Possible theoretical explanations for local maxima 
in ^(r) 

Possible theoretical explanations for fine features in the 
correlation function or power spectrum include: 

(i) acoustic baryonic fluctuations formed during the tran- 
sition of the recombination epoch, e.g. Eisenstein 

(1998) ; Meiksin, White & Peacock (1998); Peebles 

(1999) ; Miller, Nichol & Batuski (1998); 

(ii) baryonic fluctuations formed during an inflationary 
epoch, as a result of the evolution of a complex scalar 
field condensate, Kirilova & Chizhov (2000); 

(iii) phase transitions in supersymmetric models of double 
inflation, (e.g. Lesgourgues et al. 1998, 2000); or 

(iv) string theory inspired fluctuations formed during 
the Planck epoch. Brandenberger & Martin (2000); 
Easther et al. (2001) 

It would clearly be premature to attempt to distinguish 
between these on the basis of the 2QZ-10K data. 

5. Conclusions 

Although considerable care is required to avoid selec- 
tion effects in the "lOK" initial release of the 2dF QSO 
Redshift Survey (2QZ-10K), it is difficult to avoid the con- 
clusion that a harmonic of the i^gg w 130 ± 10 h^^ Mpc 
local maximum in the correlation function or power spec- 
trum of density perturbations, as traced by extragalactic 
objects, is present in the data and that this scale provides 
an extremely model-free way of constraining local cosmo- 
logical parameters without requiring combination with ex- 
ternal data sets. 

This local maximum was found by estimating the spa- 
tial two-point autocorrelation functions ^(r) of the three- 
dimensional (comoving, spatial) distribution of the N = 
2378 quasars in the most completely observed (> 85% 
spectroscopic completeness) and "covered" (> 80% cover- 
age completeness) sky regions of the catalogue, over the 
redshift ranges 0.6 < z < 1.1 ("low-z"), 1.1 < z < 1.6 
("med-2;"), and 1.6 < ^ < 2.2 ("hi-^"), using the z- 
scrambling technique in order to avoid selection effects. 



(i) Avoiding a priori estimates of the length scales of fea- 
tures, local maxima in ^(r) are found in the different 
redshift ranges. The requirement that a local maxi- 
mum be present in all three redshift ranges at the same 
comoving length scales implies strong, purely geomet- 
ric constraints on the local cosmological parameters. 
The only local maximum satisfying this requirement is 
that at a length scale of 2Llss = 244 ± 17 Mpc. 

(ii) For a standard cosmological constant FLRW model, 
the matter density and cosmological constant are con- 
strained to nm = 0.25 ± 0.10, r^A = 0.65±0.25 (68% 
confldence), Qra = 0.25 ± 0.15, f^A = 0.60±0.35 (95% 
confidence), respectively, from the 2QZ-10K alone. 

(iii) For an effective quintessence (wq) model and zero 
curvature, < —0.5 (68% confldence), wq < —0.35 
(95% confidence) are found, again from the 2QZ-10K 
alone. 

These results arc consistent with type la supcrnovae 
and microwave background results (e.g. Perlmutter et al. 
1999; Riess et al. 1998; Lange et al. 2000; Balbi et al. 2000), 
but avoid the need to combine two data sets in order to 
obtain constraints on both f2m and Qa- 

The constraints on Qm are, of course, in remarkable 
agreement with the constraints from kinematics of galaxy 
clusters (e.g. Carlberg, Yee & EUingson 1997), collapsing 
galaxy groups (Mamon 1993), as well as from the baryonic 
fraction in clusters (White et al. 1993; Mohr, Mathiesen 
& Evrard 1999) and groups (Henriksen & Mamon 1994). 

However, in contrast to Jaffe et al. (2000), who exclude 
a fiat imiverse at about ^ 95% significance, i.e. fini + 
Qa = l-lliai2 at 95% confidence, a hyperbolic universe 
is weakly favoured here, though a flat universe is consistent 
with the data at the 68% confidence level. 

The present results also contrast with the constraint 
from the quietness of the Hubble flow (Chemin 2000 and 
Sect. 5 of Sandagc 1999) which requires a cosmological 
constant of about il\ = 0.8 or higher. 

It should be emphasised that the 2QZ-10K data es- 
tablish the existence of a non-zero cosmological constant 
independently of the supernovae results: f^A = is refuted 
at the 99.7% confldence level. 

The full 2QZ will clearly provide even more impressive 
constraints on local maxima, fim, ^^A and wq. 
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